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Abstract 

Human coronavirus OC43 (HCoV-OC43) is a respiratory virus that usually causes a common cold. However, it has the 
potential to cause severe infection in young children and immunocompromised adults. Both SARS-CoV and MERS-CoV 
were shown to express proteins with the potential to evade early innate immune responses. However, the ability of HCoV- 
OC43 to antagonise the intracellular antiviral defences has not yet been investigated. The potential role of the HCoV-OC43 
structural (M and N) and accessory proteins (ns2a and ns5a) in the alteration of antiviral gene expression was investigated in 
this study. HCoV-OC43M, N, ns2a and ns5a proteins were expressed in human embryonic kidney 293 (HEK-293) cells before 
challenge with Sendai virus. The Human Antiviral Response PCR array was used to profile the antiviral gene expression in 
HEK-293 cells. Over 30 genes were downregulated in the presence of one of the HCoV-OC43 proteins, e.g. genes represent- 
ing mitogen-activated protein kinases, toll-like receptors, interferons, interleukins, and signaling transduction proteins. Our 
findings suggest that similarly to SARS-CoV and MERS-CoV, HCoV-OC43 has the ability to downregulate the transcription 
of genes critical for the activation of different antiviral signaling pathways. Further studies are needed to confirm the role of 


HCoV-OC43 structural and accessory proteins in antagonising antiviral gene expression. 


Introduction 


Human coronavirus OC43 (HCoV-OC43) is taxonomi- 
cally classified within the Coronaviridae family which also 
includes severe acute respiratory syndrome (SARS) and 
Middle East respiratory syndrome (MERS) coronaviruses. 
HCoV-OC43 belongs to subgroup A of the betacoronavi- 
ruses. These are positive-sense single-stranded RNA viruses 
with a large genome size at 30,780 bases [1]. The viral pro- 
teins that contribute to the structure of coronaviruses include 
the Membrane (M), Nucleocapsid (N), Envelope (E) and 
Spike (S) glycoprotein. Betacoronaviruses possess a unique 
protein on the surface alongside the S glycoprotein called 
Hemagglutinin esterase (HE). Each coronavirus possesses its 
own Set of unique accessory proteins. The accessory proteins 
of HCoV-OC43 are called ns2a and nsSa [2]. 
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HCoV-OC43 and other common circulating coronavi- 
ruses cause 10 to 30% of all common colds, and 5 to 10% 
of all upper and lower respiratory tract infections [3-6]. 
Subclinical or very mild infections are common, and can 
occur throughout the year. However, human coronaviruses 
other than SARS-CoV and MERS-CoV are responsible for 
10 to 20% of hospitalisations of young children and immu- 
nocompromised adults with respiratory tract illness [7]. 
HCoV-OC43 seems to have also the potential for both neu- 
rotropism and neuroinvasion that are associated with mul- 
tiple sclerosis [8, 9]. Moreover, HCoV-OC43 was detected 
in the cerebrospinal fluid of a child with acute disseminated 
encephalomyelitis [10]. A case of fatal encephalitis associ- 
ated with HCoV-OC43 infection was recently reported in an 
immunodeficient 11-month-old boy [11]. 

Several studies on human coronavirus interaction with 
host cells and the innate immune response have been con- 
ducted on SARS-CoV and MERS-CoV, because of the high 
case-fatality rate associated with these infections. In Vero E6 
cells expressing SARS-CoV N protein, interleukin 6 (IL6) 
was activated through nuclear factor kappa B (NF-«B) [12]. 
On the contrary, the M protein of SARS-CoV suppresses 
NF-«B activity and contributes to the pathogenesis of the 
virus [13]. In addition, many proteins in SARS-CoV were 
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shown to have the ability to antagonise interferon (IFN) and 
counteract the innate immune response. The ORF3b and 
ORF6 accessory proteins have the potential to antagonise 
IFN by different mechanisms [14]. SARS-CoV M protein 
was recently shown to suppress IFN production [15]. The 
structural and accessory proteins M, ORF4a, ORF4b, and 
ORF5 of MERS-CoV were shown to be potent IFN antago- 
nists [16, 17]. Itis still unknown whether the evasion mecha- 
nisms described above have a certain role in the severity 
of SARS and MERS. It is noteworthy that several people 
infected with SARS-CoV or MERS-CoV do not develop 
severe symptoms [18-21]. 

Most research on SARS and MERS has been done on 
one specific virus and how it activates or inhibits the innate 
immune response. None of the studies conducted so far have 
compared results with those of HCoV-OC43 that is asso- 
ciated with low mortality rate. HCoV-OC43 may have the 
same ability as SARS-CoV and MERS-CoV to evade innate 
immunity. Indeed, there are few studies investigating HCoV- 
OC43 interaction with host cells and the innate immune 
response. The N protein was shown to cause prolonged 
activation of NF-«B in N-expressing HEK-293T cells [22]. 
Furthermore, HCoV-OC43 infection of cell lines derived 
from human respiratory tract epithelia and hepatocytes, was 
enhanced by all three types of IFN, due to the interaction 
with the IFN-inducible transmembrane (IFITM) proteins 
[23]. Moreover, the acetyl-esterase activity of HE protein 
strongly enhances the production of infectious virus [24]. 

While SARS-CoV and MERS-CoV have received much 
attention from the scientific community, very little research 
has addressed the effect of other human coronavirus infec- 
tions on antiviral gene expression. The potential role of the 
HCoV-OC43 structural and accessory proteins in the altera- 
tion of antiviral gene expression in HEK-293 cells was inves- 
tigated in this study. 


Materials and methods 
Cell culture 


Human embryonic kidney 293 (HEK-293) cells were 
obtained from American Type Culture Collection (ATCC, 
Manassas, VA, USA), and were grown in Dulbecco’s mini- 
mal essential medium (DMEM) containing GlutaMAX™ 
(Life Technologies Corporation™, Grand Island, NY, 
USA) using 25 cm’ tissue culture flasks. The medium 
was supplemented with 10% fetal bovine serum (FBS) 
(Life Technologies ™), Fungizone® (250 ug/ml) (Life 
Technologies!™), penicillin G (10,000 U/ml) (Life Tech- 
nologies) and streptomycin sulfate (10,000 pg/ml) (Life 
Technologies‘). Monolayers of HEK-293 cell culture flasks 
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were incubated at 37 °C in the presence of 5% carbon diox- 
ide (CO,) and 90% humidity. 


Viruses 


Human coronavirus OC43 (HCoV-OC43; VR-1558), influ- 
enza A virus (H3N2, A/Hong Kong/8/68 strain; VR-544), 
and Sendai virus (VR-105) were obtained from ATCC. 


Amplification of HCoV-OC43 ns2a, ns5a, M, and N 
genes by RT-PCR 


HCoV-OC43 RNA was extracted using the QlAamp® 
Viral RNA Minikit (Qiagen™ GmbH, Hilden, Germany) 
according to the manufacturer’s instructions. The HCoV- 
OC43 ns2a, nsSa, M, and N genes were amplified by a 
two-step reverse transcription-polymerase chain reaction 
(RT-PCR) using the GcneAmp® RNA Core Kit (Applied 
Biosystems'™, Foster City, CA, USA) and 10 pmol of previ- 
ously described primers [2] on a GeneAmp® PCR System 
9700 (Applied Biosystems!™). The RT reaction conditions 
were: annealing at 37 °C for 60 minutes, denaturation at 
90 °C for 5 minutes, and cooling at 4 °C. Thereafter, the 
cDNA was amplified by PCR using the following conditions: 
denaturation at 94 °C for 10 minutes, then 35 cycles of dena- 
turation at 95 °C for 30 seconds, annealing at 60 °C for 30 
seconds, and extension at 72 °C for 30 seconds, followed by 
a final extension step at 72 °C for 7 minutes, with cooling at 
4 °C. PCR products were run on agarose gels and the bands 
that were size-specific for the amplified gene of interest were 
cut and purified by Wizard® SV Gel and PCR Clean Up 
System (Promega’™, Madison, WI, USA) according to the 
manufacturer’s instructions. Influenza A virus RNA was also 
isolated and was used to amplify the NS/ gene by RT-PCR 
using NS/-specific primers [25]. Influenza A NS1 protein 
was used as positive control for IFN antagonism [14]. 


Expression of HCoV-OC43 ns5a, ns2a, M and N 
proteins in HEK-293 cells 


The pAcGFP1-N In-Fusion® Ready Vector (Clontech™, 
Takara Bio Company, Mountain View, CA, USA) was used 
for cloning RT-PCR products. The In-Fusion® HD Cloning 
Kit (Clontech™) was used to set up the ligation reaction. 
Competent TOP10 E. coli cells (Invitrogen!™, Carlsbad, CA, 
USA) were used for transformation. The isolation of vector 
was achieved by using the PureYield'™ Plasmid Miniprep 
Kit (Promega) according to the manufacturer’s instructions. 
Confirmation of successful cloning was achieved using 
restriction digestion and direct sequencing. HEK-293 cells 
were seeded at 5 x 10° per well of 24-well plates, and trans- 
fected with ns2a-pAcGFP1, nsSa-pAcGFP1, M-pAcGFPI1, 
N-pAcGFP1, or NS1-pAcGFP1 using Lipofectamine® 2000 


Antiviral gene expression profile during HCoV-OC43 infection 


according to the manufacturer’s instructions. Expression of 
ns2a, nsSa, M, N and NS1 proteins was confirmed by indirect 
immunofluorescence assay (IFA) using anti-GFP monoclo- 
nal antibody (10 pg/ml) (Clone JL-8, Clontech™). Expres- 
sion of N and NS1 proteins was further confirmed by IFA 
using monoclonal antibody against HCoV-OC43N (EMD 
Millipore'™, Billerica, MA, USA) and polyclonal anti- 
body against influenza A NS1 protein (EMD Millipore™), 
respectively. 


Challenge of HEK-293 cells with Sendai virus 


Transfected and mock-transfected HEK-293 cells were 
grown to confluency in a 96-well plate at 37 °C in the pres- 
ence of 5% CO,. The media was then removed, and Sendai 
virus suspension (ATCC) was added to each well at 1 multi- 
plicity of infection (m.o.1). The plates were further incubated 
at 37°C with 5% CO, and at 90% humidity. 


Antiviral gene expression profile in cells expressing 
HCoV-OC43 structural and accessory proteins 


The Human Antiviral Response PCR array (Qiagen!) was 
used to profile the antiviral gene expression in HEK-293 
cells in two separate experiments. The PCR array moni- 
tored the transcriptional activity of different genes involved 
in the activation of antiviral and pro-inflammatory proteins 
following amplification of RNA transcripts by real-time 
RT-PCR. Total RNA from transfected and mock-trans- 
fected HEK-293 cells was extracted using an RNeasy® Kit 
(Qiagen™) with on-column DNase digestion according to 
the manufacturer’s instructions. The RT? First Strand Kit 
(Qiagen'™) was used for cDNA synthesis. The RT reaction 
was then added to RT? SYBR Green ROX qPCR Mastermix 
(Qiagen’™), and the PCR was carried out in a 96-well plate 
(Qiagen™) on an ABI 7500 Fast Real-Time PCR system 
(Applied Biosystems). Cycle threshold (C;) values were 
exported to an Excel file to create a table of Cy values. This 
table was then uploaded on to the Qiagen data analysis web 
portal (http://www.qiagen.com/geneglobe). Quality control 
was performed as the array contains controls to monitor 
genomic DNA contamination (GDC), first strand synthe- 
sis (RTC) and real-time PCR efficiency (PPC). Samples 
were assigned to controls and test groups. Cy values were 
normalised based on a manual selection of reference genes 
with glyceraldehyde 3-phosphate dehydrogenase (GAPDH) 
chosen as the housekeeping reference gene. Fold regulation 
comparisons were used to profile antiviral gene expression. 
Fold change/regulation was calculated using delta delta Cy 
method (AAC,) where AC; is calculated between the gene 
of interest (GOI) and GAPDH, followed by AAC, calcula- 
tions (AC, (Test Group)-AC; (Control Group)). Fold change 
is then calculated using 2“ (-AAC,) formula. Thereafter, the 


data analysis web portal was used to provide the scatter plot, 
the volcano plot, the clustergram, and the heat map. The fold 
regulation threshold was > 2 for upregulation and < -2 for 
downregulation. 


Statistical analysis 


Fold change in gene expression from 2 different experi- 
ments was summarised as mean + standard deviation. The 
difference in mean fold change between two groups was 
determined by independent-samples f-test. P-value < 0.05 
was considered significant. All statistical analyses were per- 
formed using IBM SPSS statistics software version 25.0 for 
windows (IBM Corporation, Armonk, NY, USA). 


Results 


PCR array profiling of antiviral genes in HEK-293 
cells expressing HCoV-OC43 structural 
and accessory proteins 


Figure | shows the expression of HCoV-OC43 and influenza 
proteins in HEK-293 cells. The antiviral gene expression 
in transfected HEK-293 cells expressing one HCoV-OC43 


Fig. 1 Expression of HCoV-OC43 and influenza A proteins in HEK- 
293 cells. HEK-293 cells were transfected with the expression vector 
pAcGFP1 coding for HCoV-OC43 ns2a, ns5a, M or N proteins, or 
for influenza A NS1I protein. Three days following transfection, the 
expression of viral proteins was confirmed by indirect immunofluo- 
rescence assay using anti-GFP monoclonal antibody 
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protein was compared to that in mock- and non-transfected 
cells. Only the expression of the oncogene FOS was upregu- 
lated across all samples (Table 1). Nevertheless, over 30 
genes were downregulated in the presence of one of the 
HCoV-OC43 proteins, e.g. mitogen-activated protein kinases 
(MAP3K7), toll-like receptors (TLR3 and TLRS), caspases 
(CASP/), cathepsins (CTSS), interferons (JFNA/), interleu- 
kins UL/B, ILI2A, IL15, and IL1&), and signaling trans- 
duction proteins (MYDS&& and MAVS) (Table 1). Similar 
results were obtained in the presence of the influenza A NS1 
protein, a well-known IFN antagonist [14] (Table 1). The 
mean negative fold change for CARD9 in the presence of the 
structural protein M or N was larger than that in the pres- 
ence of the accessory protein ns2a (p< 0.01). In contrast, the 
expression of the MAP3K7 gene was more downregulated 
in the presence of ns2a or ns5a protein than in the presence 
of M or N protein (p < 0.01). Downregulation of CTSS gene 
expression was greater in the presence of ns2a, ns5a or N 
protein than in the presence of M protein (p< 0.01). JL/B, 
CASP/, and TLR3 gene expression was more downregulated 


Table 1 Mean fold change 


— ns2a 

in antiviral gene expression 

in HEK-293 cells expressing Mean +SD 

HCoV-OC43 structural or 

accessory proteins FOS +9.94 0.34 
MAP3K7 -83.63 Z17 
TLR3 -55.56 0.79 
CTSS -21.91 1.64 
CASP 1 -16.73 1.56 
MEFV —— —— 
ILI2A 7.18 1.34 
PYCARD -10.06 1.24 
ILIB -24.83 0.67 
MX! -6.97 0.57 
IFIH1 -7.89 0.23 
RIPKI -9.61 1.25 
IFNAI -4.99 2.78 
CD80 -4.41 1.23 
CD86 -4.15 1.20 
CXCLI1 -5.84 0.07 
MYDS88 -7.55 ppp) 
FADD -4.61 1.23 
MAVS -3.81 1.49 
TRIM25 -4.59 1.93 
IL1I5 -3.77 0.72 
CYLD -4.98 0.04 
TLRS -2.39 0.16 
TNF -2.22 0.03 
IL18 -3.16 0.06 
TRAF6 -4.71 0.45 
IRF3 —. — 
TBK1 -3.35 0.02 
IFNARI -3.12 0.05 
SUGTI — — 
CARD9 -4.81 0.45 
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in the presence of ns2a or ns5a protein than in the presence 
of M or N protein (p < 0.01). On the contrary, CDS6 was 
more downregulated in the presence of M protein than in the 
presence of N, ns2a or ns5Sa protein (p < 0.01). 

Sendai virus, a well-known inducer of antiviral genes, 
was used to challenge HEK-293 cells. A total of 22 genes 
were upregulated in the presence of Sendai virus, whereas 
9 genes were downregulated (Table 2). There was strong 
positive fold change for MAP3K7, FOS, CXCL9, JUN, 
TLR7, IRF7, ILI2B, NOD2 and IL6 in HEK-293 cells. In 
the presence of ns2a, nsSa, M, N or NS1 protein, the chal- 
lenge of cells with Sendai virus resulted in the upregulation 
of only 7 genes and downregulation of 25 genes (Table 3). 
The expression of OAS2 was upregulated in the presence of 
HCoV-OC43 nsSa or M protein, or influenza A NS1 protein. 
A positive fold change was observed for three chemokines 
(CCL5, CXCLS, CXCL10), two RIG-I-like receptors UFTH1 
and DDX58), one interferon (JFNB1/), and one ubiquitin-like 
protein (JSG/5), whereas there was a drastic downregulation 
of MAP3K7 expression with negative fold change ranging 


nsSa M N NS1 (FluA) 
Mean +SD Mean +SD Mean +SD 

+9.49 0.23 +3.45 0.45 +10.06 0.29 +8.23 
-88.75 0 — -59.54 1.23 -83.27 
-32.59 p23 -14.68 0.63 -31.73 0.83 -25.79 
-23.85 1.85 -7.68 0.93 -24.57 0.73 -17.94 
-16.52 Los. == — -12.03 0.53 -21.32 
-14.36 152 5.21 222 21359 222 9.89 
-13.62 0.95 -6.74 1.23 -8.61 1.23 -9.59 
-12.64 0.97 -18.10 0.63 -11.48 2.56 -2.60 
-9.94 025 — — -6.63 159 -2.15 
-7.24 2 =— — -8.14 0.03 -9.12 
-6.99 1.24 -2.84 0.09 -6.05 0.94 = -7.10 
-6.90 0.75 -5.21 253 -8.85 0.35 -6.66 
-6.59 0.53 -3.24 124 -5.39 0.48 -4.42 
-6.38 0.54 -5.71 0.56  -6.04 054 — 
-6.01 0.12 -13.61 0.83 -5.69 2.34 = -4.14 
-5.05 1.26 86 -2.17 0.11 — — -2.17 
-4.79 189 -2.34 0.05 -5.15 0.53 -5.18 
-4.30 0.03 — — -5.66 0.23 -4.58 
-2.80 0.74 — — 3.15 0.05 -4.67 
-4.06 0.08 = — — -5.20 1.29 = -4.57 
-3.96 0.19 — — “3.27 0.45 -4.14 
-3.76 | Bs: se — — — -3.36 
-3.46 1.23 -3.06 0.05 -3.28 1.24 8 -2.39 
— — -7.27 0.68 = — — — 
-3.17 0.04 -3.14 0.91 -3.25 0.53 -3.14 
-3.09 0.63 -2.47 0.33 -3.76 0.73 -3.50 
-2.98 0.35 -2.19 0.03 — — -3.15 
-2.97 0.69 = -2.85 0.45 -3.05 0.23 -3.20 
-2.93 0.04 — — =2.93 0.04 -2.93 
-2.92 0.09 — — -3.97 0.05 -2.99 
— — -16.75 0.93 -12.19 0.02 — 


Results are shown as the mean + standard deviation (SD) of two independent experiments. Mean+SD was 
not available for influenza NS1 because only one experiment was performed. (——) indicates no significant 
fold change in one or two experiments 
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Antiviral gene expression profile during HCoV-OC43 infection 


Table 2 Mean fold change in antiviral gene expression in HEK-293 
cells challenged with Sendai virus 


Upregulated genes Cells (w/ Sendai) Downregu- Cells (w/ 
lated genes Sendai) 

Mean +SD Mean +SD 

MAP3K7 +276.05 8.27 JIFIHI -16.21 3.45 

FOS +177.10 2.13 CCL5 -14.95 3.43 

CXCL9 +25.64 2.71 CTSS -13.06 1.54 

JUN +2461 2.34 ILI2A -8.27 1.73 

TLR7 +21.87 1.87  TLR3 -6.53 1.35 

IRF7 +15.61 2.30 DDX58& -5.18 2.69 

IL12B +15.54 0.89 RIPKI -5.03 2.27 

NOD2 +15.26 1.74 CASPI -4.01 0.37 

IL6 +15.17 4.78 MEFV -2.29 0.01 

SPP1 +4.20 0.01 

NLRP3 +4.16 0.01 

IFNA2 $3.25 0.62 

TRADD +3.04 0.45 

MAVS +3.02 0.75 

IFNAI +2.90 0.06 

CCL3 +2.81 0.01 

IRAK1 42.78 0.12 

CTSB +2.72 0.06 

PINI 42.50 0.05 

TLRS 42.48 0.32 

IRF5 +2.36 0.06 

MAP2K3 42.31 0.21 


Results are shown as the mean+standard deviation (SD) of two inde- 
pendent experiments 


from -484 to -27,496 (Table 3). FOS, CXCL9, JUN, TLR7, 
IRF7, IL12B and IL6 also showed substantial negative fold 
changes (Table 3). The expression of MAP3K7, TLR7, 
ILI2B, IL6 and NOD2 genes was more downregulated in 
the presence of ns5a protein than in the presence of other 
HCoV-OC43 proteins (p < 0.05), whereas the expression 
of the CXCL9Y gene was more downregulated in the pres- 
ence of ns2a protein than in the presence of other proteins 
(p <0.001). Moreover, the expression of JUN and [RF7 
genes was most downregulated in the presence of N protein 
(p <0.01). Interestingly, the expression of JRF'7 and IFNA/ 
genes was more downregulated in the presence of ns2a, nsSa 
or N protein than in the presence of influenza NS1 protein. 
The expression of NOD2 was more downregulated in the 
presence of ns2a or ns5a protein than in the presence of 
influenza NS1 protein, while the expression of PYCARD was 
more downregulated in the presence of one of the tested 
HCoV-OC43 proteins than in the presence of influenza NS1 
protein. Moreover, FOS gene expression was more down- 
regulated in the presence of ns5a, M or N protein than in the 
presence of influenza NS1 protein (Table 3). 


Discussion 


The current study shows for the first time that different 
genes involved in the activation of antiviral and inflamma- 
tory responses, e.g. those representing interleukins, IFNs, 
MAPKs, PRRs and adaptor proteins, are downregulated 
in the presence of HCoV-OC43 structural or accessory 
proteins. Similar antiviral gene expression profiles were 
observed between influenza NS1 protein, which is known 
to antagonise type I IFN activity [14, 26-28], and HCoV- 
OC43 proteins. Interestingly, some genes such as [FNA/ 
and JRF7 were more downregulated in the presence of 
HCoV-OC43 proteins than in the presence of influenza 
NS1 protein. 

Accessory proteins are often not essential for virus rep- 
lication in vitro, but required for optimal replication and 
virulence in the natural host [29]. HCoV-OC43 nsS5a acces- 
sory protein is known to be a viroporin involved in virion 
morphogenesis and pathogenesis [30]. In its absence, there 
is a reduction of virulence, inflammation and replication in 
infected mice [30]. Murine coronavirus accessory protein 
5a is homologous to ns5a protein, and it has been found to 
antagonise IFN induction [31]. Our results showed that the 
expression of PRRs, JFNA/, IRF7, MYDS8, and MAVS, all 
important for the induction of an antiviral response, was 
downregulated in the presence of HCoV-OC43 accessory 
proteins ns2a and nsSa. Moreover, the expression of inter- 
leukins, chemokines, FOS, JUN, IRFs, and IFNA/ was 
downregulated in the presence of ns5a protein. FOS and 
JUN are transcription factors involved in the induction of 
inflammatory cytokines [32, 33]. Altogether, our results 
suggest HCoV-OC43 ns2a and nsS5a have the potential to 
induce a shift towards an anti-inflammatory response and 
blockage of the type I IFN pathway. However, in vivo stud- 
les are required to confirm the type of immune response 
HCoV-OC43 elicits during infection. In the absence of 
adequate immune responses, the virus will be able to rep- 
licate and spread freely. 

The M protein is indispensable in virion assembly 
between the Golgi apparatus and endoplasmic reticulum, 
and is the most abundant protein among the coronavi- 
rus proteins [34]. SARS-CoV M protein suppresses the 
activation of NF-«B [13]. It can also halt IFN production 
by blocking the formation of the TRAF3 complex [15]. 
MERS-CoV M protein is also a potent IFN antagonist, 
targeting the TRAF3-TBK1 interaction to inhibit phos- 
phorylation of IRF3 [17, 35]. Moreover, TRAF6, TRADD 
and FADD are key players in the NF-«B and apoptosis 
pathways [36]. Our results showed that the expression 
of TRAF6 was downregulated in the presence of M or 
other tested HCoV-OC43 proteins, but this downregula- 
tion could not be confirmed after stimulation with Sendai 
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Table 3 Mean fold change 


= ns2a 

in antiviral gene expression 

in Sendai-infected HEK-293 Mean +SD 

cells expressing HCoV-OC43 

structural or accessory proteins CCL5 +53.45 0.28 
IFIH1 +16.78 0.03 
IFNB1 +17.76 0.77 
ISG15 +8.68 L227 
DDX58 +14.00 0.25 
CXCL1IO —_+10.20 0.39 
OAS2 — — 
CXCL8 +9.80 1.29 
MAP3K7_ -12141.84 0.30 
FOS -12.68 1.25 
CXCL9 -104.94 0.91 
JUN -33.49 222 
TLR7 -16.23 0.01 
IRF7 -59.17 0.03 
IL12B -20.00 0.01 
NOD2 -14.32 0.01 
IL6 -12.75 0.78 
ILIB -13.65 0.33 
PYCARD ~ -6.51 0.56 
MYD88 -3.72 0.07 
CARD9 -3.71 0.01 
IKBKB -3.12 0.19 
FADD -4.87 122 
TRADD -12.70 1.05 
MAVS -9,74 0.07 
IFNA1 -18.51 0.04 
IRAKI -5.31 1.49 
CTSB —— — 
PINI -5.73 1.26 
TLR8 -3.54 0.01 
IRF5 -4.87 0.76 
MAP2K3 _ -5.34 0.10 
TICAMIL — — 


nsSa M N NS1 (FluA) 
Mean +SD Mean +SD Mean +SD 

+23.24 0.27 +30.00 0.22 +50.47 0.07 +19.74 
+8.15 0.59 415.23 0.16 +13.60 0.36 +7.47 
+7.39 0.04 +9.16 0.14 +22.13 0.22 +2.96 
+5.84 0.51 +6.87 0.70 +11.68 0.17 +3.85 
+4.80 0.59 +8.20 0.45 +10.98 0.42 +4.07 
+4.40 0.10 +8.03 0.31 +10.24 0.08 +2.91 
+3.23 0.32 +4.67 0.63 — ——- $4.52 
+3.04 0.51 +4.68 0.01 +8.88 0.71 +3.61 
-27496.30 0.44 -484.06 7.02 -16781.69 1.26 -31317.90 
-22.9] 0.60  -17.63 0.91 -29.61 2.67 -14.76 
-54.81 0.62 -46.57 3.77 -25.86 0.63 -95.64 
-51.34 1.18  -35.35 1.94 -74.12 3.17  -58.19 
-33.86 0.01 -19.42 3.91 -14.61 0.01 -53.61 
-27.66 0.72 -18.05 1.63 -87.48 3.49 = -20.75 
-41.70 0.01 -9.86 160 -13.24 0.63 -66.03 
-18.64 0.42 -5.00 2.60 -5.40 1.78  -7.68 
-37.00 0.15  -15.16 2.13 -10.19 0.51  -33.09 
-30.81 1.16 -11.86 0.16 -35.01 2.34  -28.98 
-9.06 2.71 -4.77 0.34 -7.36 1.43 -3.62 
-4.96 0.01 -3.33 0.37 — — -412 
-3.38 0.94 -3.85 160 -3.34 0.01 -3.21 
-5.14 O53 -3:57 0.07 = -7.50 2.91 -4.82 
-4.80 149 — —  £é-3.44 1.10 -5.67 
-7.15 0.94 -5.20 232 =3.91 0.83 -10.01 
-10.70 1.07 -6.03 1.13 -9.98 1.71  -10.45 
-24.28 0.20 -11.37 2.32 -21.84 1.85 -12.83 
-3.92 0.96  -5.30 1.13 -8.89 3.91 -2.47 
-7.12 0.30 -4.16 0.16 -4.91 0.17 -4.17 
-4,34 0.19 — —_  ~— — -4.94 
-7,38 0.01 -6.31 0.49 -3.18 0.01 -11.68 
-4.14 0.11 -5.32 115  -5.86 2.59 -2.98 
-10.28 0.27 -8.57 0.10 -4.83 0.42 -8.88 
-3.51 0.69 -7.15 1.74 — — ~=«=-3.51 
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Results are shown as mean+standard deviation (SD) of two independent experiments. Mean+SD was 
not available for influenza NS1 as only one experiment was performed. (——) indicates no significant fold 


change in one or two experiments 


virus. In contrast, the downregulation of TRADD was only 
observed after challenge with Sendai virus. However, the 
expression of FADD was not downregulated in the pres- 
ence of HCoV-OC43M protein before or after challenge 
with Sendai virus. Furthermore, the expression of MAVS, 
MYDS&%&8, IRF and IFNA/ that are essential in the induction 
of the type I IFN signaling pathway, was downregulated 
in the presence of M protein. This suggests that HCoV- 
OC43M protein has similar downregulatory effects on the 
induction of type I IFN and NF-«B as SARS-CoV and 
MERS-CoV M protein. Further studies are needed to con- 
firm the effects of HCoV-OC43M protein on the type I IFN 
and NF-«B signaling pathways. 

The function of the N protein is to ensure the replica- 
tion and transcription of viral RNA while disrupting the cell 
cycle [37, 38]. Indeed, SARS-CoV N protein inhibits the 
progression of the cell cycle and activates the proinflamma- 
tory factor cyclooxygenase-2 [39]. Moreover, it was shown 
that SARS-CoV N protein inhibits JFN-f expression [14], 
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while it activates JL-6 expression [13]. HCoV-OC43N pro- 
tein potentiates NF-«B by binding its inhibitor, microRNA 
9 [22]. Our study did not show significant fold changes in 
NF-«B expression in the presence of N protein; however 
there was a significant downregulation in expression of 
PRRs (NOD2 and TLR7), adaptor proteins (MYD&8, MAVS, 
CARD9Y) and signal mediators UKBKB, IRAKI, TBK1, 
RIPKI, TRADD, TRAF), all critical in the induction of 
NF-«B signaling pathway. Moreover, our results showed 
that expression of MAPKs (MAP3K7 and MAP2K3) was 
downregulated in the presence of N protein, which is in line 
with a previous study showing that murine coronavirus N 
protein inhibits signaling of the AP-1 complex [40], which 
is composed of FOS and JUN oncogenes, whose expression 
was also downregulated in our investigation - in the presence 
of HCoV-OC43 ns2a, ns5a, N or M proteins after challenge 
with Sendai virus. 

MAP3K7, also called Transforming growth factor-b (TGF- 
b)-activated kinase I (TAK1), was the most downregulated 
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gene in the presence of any HCoV-OC43 protein tested 
after challenging cells with Sendai virus. TAK] 1s a protein 
kinase known to be a critical mediator of the inflammatory 
response as it is activated by lipopolysaccharide (LPS) and 
other TLR/NLR agonists, TNF-a and IL-1 [41]. Binding of 
IL-1 and LPS to their receptors activates TAK1 through a 
common pathway whereby the kinases IRAK1 and [IRAK4 
recruit TRAF6. TNF-a activates TAK1 via a similar mecha- 
nism, which is mediated by RIP1 and TRAF2. Activation of 
TRAF2 and TRAF6 leads to the generation of K63-linked 
polyubiquitin chains, which bind to the adaptor proteins 
TAB2 or TAB3. TAK] is bound constitutively to the acces- 
sory protein TABI and to the homologs TAB2 or TAB3. 
Once activated, TAK1 transduces the signal to NF-«B, 
c-Jun N-terminal kinase (JNK), and p38 via phosphoryla- 
tion of IKK, mitogen-activated protein kinase kinase 4/7 
(MKK4/7), and MKK3/6, respectively. Ultimately, NF-«B 
and other transcription factors downstream of p38 and JNK 
are activated, resulting in the transcription of genes impor- 
tant for inflammatory and immune responses. Interestingly, 
Protein kinase R (PKR) interacts with TRAF6, TAK1, and 
TAB2 as a scaffold protein to mediate TLR3-dependent 
NF-«B and MAPK activation [42, 43]. The strong down- 
regulation of TAK1 expression observed in our study could 
be due to: 1) the inhibition of TAK1 expression by direct 
interaction of HCoV-OC43 protein with TAK1 preventing 
binding of other mediators like TRAF2/6 or TAB2/3; 2) the 
upregulation of expression of negative regulators of TAK1 
like PP2C or PP6 phosphatases, and deubiquitinases such as 
ubiquitin-specific peptidase 4 (USP4) [41]. 

In certain cell lines, IFITM proteins were shown to 
facilitate HCoV-OC43 entry into host cells [23]. Our study 
showed that all tested HCoV-OC43 proteins were able to 
downregulate the expression of IFNA1, suggesting resist- 
ance to IFN-a induction. The interaction between HCoV- 
OC43 proteins and the host cell requires further investigation 
to understand the processes that could contribute to overall 
resistance to the antiviral state. Having all these HCoV- 
OC43 proteins with potential inhibitory effects expressed 
in close proximity in infected cells indicates a high replica- 
tive capacity of the virus. 

In conclusion, our findings suggest that similarly to 
SARS-CoV and MERS-CoV, HCoV-OC43 has the ability 
to downregulate the transcription of genes critical for the 
activation of different antiviral signaling pathways. While 
SARS-CoV and MERS-CoV were shown to induce proin- 
flammatory cytokines in vitro [44, 45], HCoV-OC43 acces- 
sory and structural proteins have the potential to inhibit 
inflammatory response in HEK-293 cells. Further studies are 
in progress to confirm the role of HCoV-OC43 structural and 
accessory proteins in antagonising antiviral gene expression, 
and to elucidate the nature of the immune response and the 
level of inflammation seen during HCoV-OC43 infection. 
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